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1.0 


INTRODUCTION 


Since  mid  1976,  the  Space  and  Missile  Test  Center  (SAMTEC)  at  Vandenberg  Air 
Force  Base  has  been  investigating  a  technique  called  Range  Vernier  whereby 
radial  range  measurements  are  derived  from  the  phase  of  signals  that  are 
received  during  operations  of  certain  pulse  Doppler  radars  on  the  Western  Test 
Range  (WTR).  These  radars,  termed  CSP  radars,  perform  fine  line  tracking  of 
Doppler  frequencies.  Implementation  of  Range  Vernier  involves  extraction  and 
recording  of  Doppler  phase  measurements,  and  subsequent  conversion  to  radial 
range  measurements.  The  latter  measurements  are  called  Doppler  ranges. 

The  theoretical  quality  of  extreme  precision  plus  modelable  error  components 
to  be  provided  by  Doppler  ranges,  offers  an  opportunity  to  greatly  improve 
accuracies  of  WTR  tracking  instrumentation.  In  addition  to  improving  the 
quality  of  range  data  per  se,  the  excellent  precision  of  the  range  measure¬ 
ments  is  especially  suited  to  substantial  improvement  of  (a)  trajectory  re¬ 
construction,  (b)  instrumentation  performance  analysis,  (c)  weapon  system 
analysis,  and  in  some  cases,  (d)  refinement  of  geodetic,  geopotential,  and 
meteorological  models.  Each  of  those  factors  is  fundamental  to  the  principal 
measurement  objectives  of  test  mission  support. 

The  purpose  of  this  document  is  to  support  further  SAMTEC  development  of  Range 
Vernier  capabilities  by  presenting  technical  information  that  is  relevant  to 
the  technique  in  general,  and  that  is  directly  pertinent  to  implementation  of 
the  technique  on  CSP  radars  at  the  WTR^.  A  major  portion  of  the  technical 
material  presented  herein  has  been  collected  from  the  documents  cited  in  the 
References,  page  34,  especially  the  paper  by  Dr.  R.  A.  Brooks,  Reference  [1], 

The  following  sections  include  a  background  and  conceptual  definition  of  the 
Range  Vernier  technique,  discussion  of  Doppler  cycle  ambiguity  resolution, 
mechanization  and  requirements  of  Doppler  phase  measurements,  analysis  of 
Doppler  range  errors,  various  applications  of  Doppler  ranges,  a  discussion  of 
feasibility  testing,  and  a  summary. 


''"A  SAMTEC  task  to  implement  a  prototype  Range  Vernier  subsystem  on  the  F PQ- 6 
CSP  radar  at  Pillar  Point,  California  is  currently  in  progress. 
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2.0 


RANGE  VERNIER  CONCEPT 


Fundamentals  of  a  certain  technique  whereby  Doppler  frequency  measurements  are 
combined  with  small  fractional  measurements  of  Doppler  phase  angles  to  subse¬ 
quently  derive  very  precise  measurements  of  radial  ranges  are  described  in  this 
section.  The  technique,  called  "Range  Vernier",  is  intended  for  application 
to  data  provided  by  a  pulse  Doppler  tracking  radar  (i.e.  a  CSP  radar)  that  is 
augmented  with  Doppler  phase  measurement  capabilities. 

2.1  Basic  Relation 

To  establish  a  basic  relation  between  the  Doppler  phase  ©d(t)  of  a  signal  reply 

to  the  radar  at  time  t  and  the  range  R  from  the  radar  to  the  object  being  tracked, 

consider  a  radar  that  is  operating  with  a  stable  transmitter  frequency  f  ,  and 

assume  there  is  no  time  delay  of  this  signal  at  the  target.  Consequently,  the 

phase  0Q  of  this  coherent  signal  is  varying  according  to  the  constant  rate  2nf  , 

and  the  difference  0Q(t  )  -  ©o(tb)  of  this  transmitter  signal  phase  between 

times  t  and  t.  is 
3  b 

VV  ■  VV  =  2nfo<-\  -  V  ■  <21> 

Also,  let  R( t ' )  be  the  range  corresponding  to  t 1 ,  and  let  t1  be  the  one-way 
transit  time,  so  that  the  time  t  of  signal  reply  at  the  radar  is  t  =  t 1  +  x1. 

The  Doppler  phase  ©d(t)  of  this  signal  reply  is  defined  by  the  equation 

0d(t)  =  0Q(t)  -  0o(t  -  2x')  . 

Thus,  because  of  equation  (2.1), 

0d(t)  =  2/tfo[t  -  (t  -  2t')]  =  47ifox'  . 

The  preceeding  relation  is  equivalent  to  the  expression 

0d(t)  =  TR(t,)  »  (2.2) 


2 


where  A  =  C/f  is  the  wavelength;  R(t')  =  Cx ‘ ;  and  C  is  the  speed  of  light 
through  the  medium. 

In  other  words  the  range  R(t')  at  the  retarded  time  t'  is  equal  to  the  scalar 
A/4 n  times  the  Doppler  phase  of  the  signal  reply  at  time  t. 

2.2  Range  Measurement  Construction 

During  operation,  a  pulse  Doppler  tracking  radar  provides  a  continuous  reference 
signal  with  constant  frequency  f  .  The  transmitted  signal  however,  consists 
of  a  train  of  pulse  modulated  bursts  of  the  reference  signal.  The  bursts  are 
separated  in  time  according  to  some  specific  pulse  repetition  frequency  f  . 
Therefore,  direct  measurement  of  Doppler  phase  between  pulses  is  not  possible 
due  to  absence  of  signal  replies  between  pulses.  However,  consideration  of 
the  true  delta-phase  D^  between  successive  pulses  is  practical;  that  is 
consideration  of  the  difference 

°i  =  W3,' 

where 

Vi  =  ed(ti+i)  and  ei  =  w  • 

A  vital  point  of  the  Range  Vernier  technique  is  that  the  signal  present  during 
a  pulse  reply  can  be  passed  through  a  measurement  subsystem  (discussed  in 
Section  4.0)  to  obtain  certain  angular  measurements  0..  These  are  small  fractions 
of  the  total  Doppler  phases  0^.  More  specifically,  each  0^  is  a  measurement 
of  the  true  Doppler  phase  0^,  modulo  2n  plus  an  error  e^.  That  is, 

0i  =  01  -  2Kin  +  e.  ,  (2.3) 

where  is  a  positive  integer  such  that  0  <  0^  2Kin  <  2n  and  is  the  error 
of  0^. 


Notice  that 


where 


ei+l  - 

0i  = 

Di  - 

Ni  = 

Ki+1- 

Ki 

(2.4) 


N.j  is  the  count  (or  multiple)  of  complete  Doppler  phase  cycles  between  successive 
pulse  replies.  The  integer  can  be  determined  with  an  ambiguity  resolution 
technique  described  later  in  Section  3.0.  Accordingly,  the  defined  below 
is  a  measurement  of  D^ 


°i  = 


0i+l  ‘  6i 


2N.ti 


(2.5) 


Also,  using  equation  (2.4),  an  equivalent  expression  is 


Di  =  °i  *  £i+i  ■  ei 


In  order  to  continue  construction  of  ranges,  form  the  sum 


S  =  D  +  D,  +  •  •  •  +  D_  , 
n  o  1  n-1 

of  delta-phase  measurements  D. ,  0  <  i  <  (n-1). 

Straight  forward  calculation  of  the  right  hand  side  of  Sn  shows  that 

S  =  0  -  0  +  e  -  e  . 

n  n  o  n  o 

A  key  point  is  that  due  to  mutual  cancellation,  the  Doppler  phases  0.  and  the 

»  J 

errors  e.  for  0  <  j  <  n  vanish  in  calculation  of  the  sum.  This  feature,  zero 

J 

accumulation  of  errors  is  very  important  to  precision  of  the  range  con¬ 
struction  with  the  Range  Vernier  technique.  To  continue  explanation  of  the 
range  construction,  the  measurement  ARr  is  defined  by  the  equality 


AR„  =  aS„ 
n  n 
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where  a  =  (\/4n)  and  \  is  the  wavelength  of  the  reference  signal.  Recall,  from 
the  basic  relation  of  equation  (2.2),  that  a  is  the  scalar  needed  for  conversion 
of  Doppler  phase  angles  to  range  units,  therefore 


AR„  =  R  -  R  +  ae„  -  ae  . 
n  n  o  n  o 

A  range  ambiguity  exists  at  this  stage  of  construction  due  to  the  unknown  value 

Rq;  but,  in  practice,  a  measurement  Rq  of  RQ  can  be  obtained  from  the  range 

subsystem  of  the  radar  (or  from  other  measurement  sources)  such  that  Rq  =  R^  + 

E  .  where  E  is  the  error  of  R  . 
oo  o 

Accordingly,  the  calculation 

R  =  ARW  +  Rrt 
n  No 

determines  a  range  R  constructed  with  the  Range  Vernier  technique.  Ranges 
determined  in  this  manner  are  called  "Doppler  ranges". 

Additional  calculations  lead  to  the  expression 

R„  =  R  +  B  +  ae  ,  where  =  E  -  as  . 
n  non  o  o  o 

This  indicates  that  as  n  increases,  the  errors  (R  •  R  )  consist  of  two  com- 

’  n  n 

ponents,  a  constant  bias  term  Bq  and  a  time  varying  error  ofen-  A  more  complete 
discussion  of  Doppler  range  errors  is  presented  in  Section  5.0. 

Measurements  from  sources  other  than  Doppler  phase  information  are  needed  to 
determine  the  multiple  Ni  of  complete  Doppler  phase  cycles  between  successive 
pulse  replies  that  is  shown  in  equation  (2.5)  above.  A  technique  that  can  be 
used  to  determine  N.,  and  consequently  resolve  Doppler  phase  ambiguity  is 
described  in  the  next  section. 
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3.0 


AMBIGUITY  RESOLUTION 


The  construction  of  range  measurements  determined  with  the  Range  Vernier 
technique  described  above  in  Section  2.2,  assumes  that  ambiguities  of  the  delta- 
phase  measurements  have  been  resolved.  These  particular  ambiguities  are  due 
to  the  unknown  numbers  of  complete  cycles  of  the  Doppler  phase  that  occur 
between  successive  signal  replies.  This  section  herein  describes  a  method  to 
resolve  delta-phase  ambiguities  by  determination  of  the  integers  such  that 


A  conceptual  background  and  general  conditions  that  are  sufficient  to  resolution 
of  delta-phase  ambiguity  are  presented  in  Section  3.1. 

A  particular  technique  to  resolve  delta-phase  ambiguity  is  described  in  Section 
3.2.  The  technique  combines  Doppler  phase  measurements  with  conventional  Doppler 
frequency  measurements.  Also,  the  technique  assumes  that  ambiguities  of  the 
latter  frequency  measurements  have  been  resolved. 

In  order  to  complete  this  discussion  of  ambiguity  resolution,  several  methods 
that  can  provide  unambiguous  Doppler  frequencies  are  mentioned  in  Section  3.3. 

3.1  Ambiguity  Background 

The  problem  of  ambiguity  resolution  can  be  addressed  by  consideration  of  the 
fundamental  derivative  expression  that  relates  Doppler  phase  ©d  to  Doppler 
frequency,  namely 


x  d  ed(t) 

2 n  ~dt 


Fd(t) 


(3.1) 
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For  the  continuous  case,  (3.1)  can  be  integrated  to  obtain  the  expression 

ti+l 

Di  =  ed(ti+l5  "  ed(tiJ  ~  Zn  J"  Fd(t)dt  .  where 

ti 

Where  is  the  delta-phase  between  t^  and  t^+,.  The  integral  on  the 
right  hand  side  can  be  expressed  as  an  integer  Ni  plus  a  non-negative 
fraction  that  is  less  than  1  to  obtain  the  expression 

Di  =  2n  (Ni  +  qi )  . 

This  shows  that  the  number  of  complete  cycles  of  Doppler  phase  between 
times  t.  and  t.+-|  can  be  determined  by  integration  of  the  Doppler  frequencies. 

However,  in  contrast  to  the  continuous  case,  operation  of  a  CSP  radar  provides 
the  Doppler  frequency  measurement  F d ( t ^ )  at  discrete  times  ti  according  to  a 
specific  pulse  repetition  frequency  f  .  As  a  consequence,  the  integral  of 
Doppler  frequencies  must  be  approximated  between  successive  pulse  replies 
that  are  separated  by  1 /f r  units  of  time;  i.e.,  t^+1  -  t.  =  1/f  .  Integral 
approximation  is  not  difficult  for  customary  test  operations  at  the  WTR 
because  in  test  missions,  Doppler  frequencies  vary  slowly  and  very  good 
approximations  to  the  integral  can  be  made,  provided  that  ambiguities  are 
not  present  in  the  measurements  F^(t. ). 

To  continue  this  background  discussion  of  the  method  described  in 
Section  3.2  below,  consider  the  following  general  case. 


Assume  the  D^  is  an  estimate  of  Di  that  has  been  previously  determined, 
for  instance  by  approximation  of  the  Doppler  frequency  integral  or  from 

some  arbitrary  source.  Also,  let  6D^  be  the  error  of  D^. 

Let  the  delta-phase  measurement  D.  be  defined  by 


(3.2) 
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where  0.^  and  0^  have  been  obtained  from  a  Doppler  phase  measurement  sub¬ 
system  and  is  an  integer  such  that 


For  this  measurement  of  D^  and  selection  of  L.,  a  sufficient  condition  for  correct 
resolution  of  the  delta-phase  ambiguity  is  the  inequality 


I  6Di  "  <ei+i  -  |<  n  .  (3.3) 

This  condition  is  not  difficult  to  meet  in  practice  whenever  Doppler  frequencies 
have  been  resolved,  because  with  operational  CSP  radars  and  f  =  160  Hz,  6Di  is 
less  than  2.25  degrees  and  (e.+^  -  e.)  is  less  than  15  degrees. 

For  this  example,  |  6D.  -  (e.+^  _  e^))  <  17.25  degrees.  This  would  certainly 
satisfy  the  condition  for  resolution  of  the  delta-phase  ambiguity.  The  above 
background  shows  the  need  of  a  solution  for  the  unknown  value  of  the  integer 
L..  A  procedure  whereby  can  be  evaluated  is  described  in  the  next  sub¬ 
section. 


3.2  Delta-Phase  Ambiguity  Resolution 


The  ambiguities  of  Doppler  frequency  measurements  F^Ct^)  are  assumed  to  have 
been  correctly  resolved  in  this  subsection.  With  this  assumption,  the  calcu- 
lational  procedure  to  resolve  delta-phase  ambiguities  can  be  described  as  follows. 

Integrate1  the  two  most  recent  Doppler  frequencies  Fd(t..)  and  Fd(t.j+1)  corres¬ 
ponding  to  successive  pulse  replies  to  obtain  the  integer  Ni  and  fractional 
value  q.j  such  that 


1 


A  trapezoidal  rule  will  provide  the  integral  approximation 

Si  .  W  -  f>m>  J  where ,  .  . ... 

2n  _  T  ’  where  ti+l  t1  ~  1/fr 
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D.  ti+l  ~  . 

2~  =  /  Fd(t)dt  =  N.  +  qi  where  0  <  q^  <  1 


Then  determine  the  comparison  value 


where  di  =  6i+1  -  6i  and  the  direct  phase  measurements  0.+1,  0i  are  obtained 
from  the  (I/Q)  subsystem. 

Notice  the  range  of  is  between  -1  and  2  since  |  d.-  |  <  2n  and  0  <q.  <  1. 

Select  the  appropriate  value  of  AN^  according  to  the  functional  list  and  range 
of  values  below 

ANi 

-1  -1  <  C.  <  -1/2 

0  -1/2  <  C.  <  1/2 

1  1/2  <  Ci  <  3/2 

2  3/2  <  Ci  <  2 

Calculate  =  N.  +  AN.,  and 

Di  =  d.  +  2N.n  ,  where  N.  =  L.  (3.4) 

As  a  consequence  of  this  procedure  the  following  inequality  is  satisfied 
7i  >  |  D.  -  D.  |  =  |  2(N.  -  L.)n  +  2 nq.  -  d  | 
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A  point  of  interest  here  is  that  delta-phase  ambiguities  can  be  resolved  for 
successive  time  instants  t. ,  i  =0,  1,  ....  N  provided  that  the  measurements 
D.  of  equation  (3.4)  can  be  evaluated  and  the  inequality  of  (3.3)  is  satisfied. 
Subsequently,  the  Doppler  ranges  Rn  corresponding  to  retarded  times  t^, 
n  =  1,  2,  ...,  N  can  be  derived  with  the  construction  of  Section  2.2.  However, 
in  practice,  there  may  be  periods  of  time  when  the  data  are  not  acceptable  or 
unavailable  for  calculation  of  Doppler  ranges;  e.g.  no  signal  replies,  gross 
errors,  lack  of  unambiguous  Doppler  frequencies,  etc.  Whenever  this  occurs, 
the  Range  Vernier  technique  must  be  reinitialized  and  started  at  a  new  time 
instant  t*  that  is  followed  by  successive  time  instants  that  do  provide 
acceptable  data.  This  implies  a  new  range  measurement  RQ(t£)  must  be  pro¬ 
vided,  and  that  another  estimate  of  the  bias  Bq,  of  Section  2. ,0  is  anticipated. 

Due  to  the  vital  role  that  unambiguous  Doppler  frequencies  play  in  the  resolu¬ 
tion  of  delta-phase  ambiguities,  a  brief  discussion  of  Doppler  frequency 
ambiguity  is  presented  in  the  next  section. 

3.3  Doppler  Frequency  Ambiguity 

Doppler  frequency  ambiguities  are  considered  to  be  resolved  whenever  the 
magnitudes  of  Doppler  frequency  errors  are  less  than  half  of  the  pulse 
repetition  frequency.  This  condition  can  be  attained  with  several  methods 
discussed  below. 

One  method  is  to  compare  (a)  range  measurements  collected  by  a  range  tracking 
subsystem  with  (b)  measurements  of  range  derived  by  integration  of  range  rates 
obtained  from  Doppler  frequency  measurements  collected  by  the  Doppler  tracking 
subsystem.  The  comparisons  can  be  processed  with  a  linear  estimator,  such  as 
a  Kalman  type  recursive  filter  or  a  batch  least  squares  filter,  to  resolve 
Doppler  frequency  ambiguitites.  A  particular  example  of  this  method  is  described 
in  the  Appendix. 

Another  method  involves  use  of  an  Invariant  Embedding  technique  to  resolve 
frequency  ambiguities.  This  method  is  currently  available  with  CSP  radars. 

It  is  described  in  Reference  [2]. 
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Either  of  the  two  above  methods  is  feasible  for  realtime  or  post-flight  data 
reduction  purpose.  However,  a  comparison  of  the  methods  indicates  that  the 
former  is  suitable  for  more  rapid  ambiguity  resolution  than  the  latter.  Rapid 
ambiguity  resolution  is  crucial  to  realtime  applications,  so  the  former  method 
is  a  strong  candidate  for  future  implementation. 

A  third  method  is  to  transform  coordinates  of  a  suitably  accurate  estimate  of 
the  trajectory  to  measurements  of  Doppler  frequencies  that  are  unambiguous. 
More  generally,  this  can  be  done  in  post-flight  applications,  when  an  accurate 
best  estimate  of  trajectory  (BET)  is  available. 

Accordingly,  there  are  several  methods  whereby  Doppler  frequency  ambiguities 
can  be  resolved  as  required  by  the  delta-phase  ambiguity  resolution  that  is 
described  in  Section  3.2. 
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4.0 


MECHANIZATION 


Mechanization  of  Doppler  phase  measurements  at  Coherent  Signal  Processor  (CSP) 
Radars  of  the  Western  Test  Range  (WTR)  is  discussed  in  this  section.  Essentially, 
extraction  and  recording  of  Doppler  phase  and  frequency  is  involved. 

4.1  Background  Discussion 

A  block  diagram  of  a  SAMTEC  Doppler  radar  (CSP  radar  of  interest)  displaying 
the  signals  and  nominal  frequencies  at  various  system  modes,  is  illustrated  in 
Figure  4.1.  The  stable  clock  and  the  frequency  synthesizer  are  external  to 
the  tracking  loop,  but  provide  the  reference  signals  required  for  tracking  loop 
operation  as  well  as  the  transmitter  carrier  signal  e^o*. 

The  pulse  amplitude  modulated  received  signal  arrives  within  the  range  gate 
window,  and  its  carrier  e^U)ot  0^  contains  the  Doppler  phase  modulation  0. 

Internal  to  the  tracking  loop,  a  digital  frequency  synthesizer  (DFS)  develops 
a  signal  e^wL^  0  ^  with  phase  modulation  0'  approximately  equal,  modulo  2n, 
to  the  Doppler  phase  0.  This  signal  is  translated  in  frequency  and  heterodyned 
with  the  received  signal  to  obtain  the  IF  input  signal.  The  IF  bandwidth  is 
approximately  matched  to  the  received  pulsewidth,  and  consequently  the  SNR  at 
the  IF  output  is  near  optimum.  The  IF  output  signal  carrier  is  ei[u»jt+A0]^ 
where  A0  =  0  -  0‘  ,  and  thus  this  signal  together  with  the  digital  frequency 
synthesizer  signal  e^U)Lt  0  ^  jointly  contain  the  Doppler  phase  information. 

The  conventional  Doppler  frequency  output  is  shown  as  F^,  but  the  remainder  of 
the  tracking  loop  signals  and  operations  are  not  germane  to  the  current  dis¬ 
cussion  and  will  not  be  considered  further. 

4.2  Doppler  Phase  Measurement 

Doppler  phase  measurement  could  be  accomplished  at  the  input  of  the  Doppler 
tracker  by  direct  comparison  of  the  phases  of  the  received  signal  and  trans¬ 
mitter  carrier,  but  there  are  several  reasons  why  this  technique  may  not  be 
preferred. 
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SYNTHESIZER 


First,  in  order  to  enhance  the  SNR,  it  is  necessary  to  filter  the  received 
signal  prior  to  measuring  its  phase.  For  near  optimum  SNR,  the  filter  band¬ 
width  should  be  matched  to  the  received  pulsewidth.  But  the  filter  bandwidth 
must  also  be  large  enough  to  accommodate  the  band  of  Doppler  frequencies,  and 
moreover,  unless  the  bandwidth  greatly  exceeds  this  band,  the  filter  will 
introduce  phase  error  in  amounts  related  to  the  instantaneous  Doppler 
frequency.  Thus,  some  compromise  in  SNR  may  result  from  a  filter  designed  to 
have  a  constant  phase  response  over  the  Doppler  frequency  band. 

Second,  when  the  Doppler  frequency  is  of  sufficient  magnitude  to  produce  a 
substantial  change  in  the  Doppler  phase  within  the  received  pulse,  precise 
timing  of  the  phase  measurement  is  required.  Thus  direct  comparison  of  the 
received  signal  and  transmitter  carrier  phases  requires  measurement  at  a 
precisely  timed  point  within  the  received  pulse.  This  can  be  accomplished 
with  high  speed  circuitry,  but  since  only  a  portion  of  the  received  pulse  is 
used,  some  information  is  lost  and  measurement  performance  is  degraded.  It  is 
also  possible  to  compute  the  average  phase  over  the  entire  pulse,  but  then 
uncertainty  in  the  time  tag  to  be  applied  to  the  measurement  is  created  by 
intrapulse  amplitude  fluctuations  as  well  as  pulse  time  of  arrival  error 
induced  by  noise. 

Third,  in  systems  which  employ  pulse  coding,  it  is  necessary  to  pass  the 
received  pulse  through  a  matched  filter  prior  to  extracting  the  Doppler  phase. 
The  performance  of  the  matched  filter  may  be  enhanced  by,  or  in  some  cases 
critically  dependent  upon,  Doppler  frequency  and  time  dilation  compensation. 
This  compensation  requires  closed  loop  operation  which  is  not  consistent  with 
the  direct  measurement  of  Doppler  phase  by  simply  comparing  the  phases  of  the 
received  signal  and  transmitter  carrier. 

An  indirect  method  of  Doppler  phase  measurement  is  conceptually  illustrated  in 
Figure  4.2.  In  the  upper  channel  the  phase  of  each  received  pulse  is  measured 
at  the  IF  output,  and  in  the  lower  channel  the  phase  of  the  DFS  output  signal 
is  also  measured.  The  phase  measurements  in  these  two  channels  are  then 
summed  to  obtain  the  Doppler  phase  measurement. 
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When  the  system  is  tracking  properly,  the  phase  variation  within  each  pulse  in 
the  upper  channel  is  negligible.  Thus  the  average  phase  over  the  pulse  is 
measured,  and  precise  timing  of  the  measurement  is  unnecessary.  On  the  other 
hand,  the  phase  of  the  DFS  signal  is  varying  at  a  rate  approximately  equal  to 
instantaneous  Doppler  angular  frequency,  and  precise  timing  of  the  phase 
measurement  in  the  lower  channel  is  required.  However,  the  DFS  output  signal 
is  continuous  and  noisefree,  and  consequently  precise  measurements  of  phase 
and  time  at  any  point  for  this  signal  can  be  obtained  with  readily  available 
hardware  components.  It  is  important  to  note  that,  although  the  time  of  the 
phase  measurement  must  be  precisely  measured,  the  actual  time  at  which  the 
phase  measurement  occurs  is  not  too  critical;  generally  it  suffices  for  the 
measurement  time  to  be  anywhere  within  or  very  near  the  range  gate. 

A  specific  mechanization  for  the  indirect  method  of  Doppler  phase  measurement 
is  illustrated  in  Figure  4.3,  and  two  options  are  presented  for  the  measure¬ 
ment  of  the  DFS  output  signal. 

The  IF  output  signal  is  passed  through  a  saturating  amplifier  and  harmonic  re¬ 
jection  filter  prior  to  phase  measurement  in  order  to  normalize  the  signal 
level  at  the  input  of  the  I/Q  device.  The  bandwidths  of  the  amplifier  and 
filter  must  be  greater  than  the  IF  bandwidth  in  order  to  not  distort  the  IF 
phase  information.  However,  the  harmonic  rejection  filter  should  provide  at 
least  40  dB  attenuation  of  all  IF  harmonics  in  order  to  ensure  negligible 
mechanization  error.  The  I/Q  ouputs  are  averaged  over  the  range  gate  window 

to  obtain  I.  and  QA. 

A  A 

In  the  first  option  for  DFS  output  processing,  a  simple  I/Q  device  is  used, 
and  its  outputs  are  sampled  at  precise  times  which  can  be  controlled  by  a 
clock  driven  trigger  within  the  range  gate. 

In  the  second  option,  a  commercially  available  computing  counter  is  used  to 
precisely  measure  the  time  interval  between  positive  slope  ze^o  crossings  of 
the  clock  signal  and  the  DFS  output  signal  (after  translation  to  the  clock 
frequency).  The  arming  signal  is  generated  by  a  negative  slope  zero  crossing 


OUTPUT  PROCESSING 


ARMING  SIGNAL 


of  the  clock  signal  selected  to  be  within  the  range  gate.  Thus  the  time  t_  of 
the  arming  signal  is  known,  and  the  computing  counter  measures  the  time  interval 
(with  proper  sign)  between  the  first  positive  slope  zero  crossing  to  follow 
arming.  Using  the  known  time  of  the  arming  command  t_  and  the  measured  time 
interval  TI,  it  follows  that  0'  =  -  u>c(TI)  at 


t  =  t  +  — 
-  10. 


+  TI 


and  thus  the  measured  Doppler  phase  at  time  t  is  simply 
0(t)  =  A0  -  w  (TI). 

The  rms  timing  precision  achievable  by  this  technique  is  better  than  0.1  nano¬ 
second. 


4.3 


Mechanization  Requirement 


A  listing  of  mechanization  requirements  pertinent  to  implementation  of  the  I/Q 
measurement  sub-system  is  presented  in  Table  4.1.  Also,  a  historical  record 
of  the  Doppler  frequency  measurements  corresponding  to  successive  replies  at 
times  t.j ,  t..+p  ...  is  required.  These  frequencies  should  be  recorded  at  the 
pulse  repetition  rate,  i.e.  f  per  second. 
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TABLE  4.1  MECHANIZATION  REQUIREMENTS 


RECORDING 

RATE 

I' ,  Q1  (OR  0') 

IA,  Qa  (OR  A0) 

PRF 

PRF 

SYNCHRONIZATION 

COINCIDENT  WITH  RANGE 

GATE  TO  NEAREST  1  pSEC 

COINCIDENT  WITH 

RANGE  GATE 

TIME  TAG 

LSB 

JITTER  (1  ct) 

1  pSEC 
<  0.05  pSEC 

N/A 

QUANTIZATION 

WORD  LENGTH* 

8  BITS 

8  BITS 

SYSTEM  MECHANIZATION 
ERROR 


FULL  WORD  REPRESENTS  EITHER  360°  IN  0  OR  1  IN  I,  Q 
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5.0 


ERROR  ANALYSIS 


A  general  model  of  the  error  present  in  range  measurements  determined  with 
the  Range  Vernier  technique  was  mentioned  in  Section  2.2.  The  general  model 
is  expanded  herein  to  recognize  more  components  of  the  basic  error  sources. 

5.1  Range  Vernier  Error  Model 

The  expanded  error  model  is 


R^  -  R.  +  6B  +  (R^-RQ)6t  +  (R.-RQ)6F  +  6p^  -  6pQ  +611.  -  6Uq  +  ,  i<0 

where 

6B  -  bias  due  to  initialization 

6t  -  systematic  time  tag  delay  error 

6F  -  scale  factor  error  due  to  uncertainties 

in  the  speed  of  light  and  long  term  clock  frequency 

6p  -  refraction  error 

6U  -  structurally  unspecified  errors  due  to  transmitter/ 

receiver  mechanization,  vehicle  transponder/antenna/plume 
and  scintillation 

v  -  serially  uncorrelated  phase  measurement  error  due  to  noise, 
jitter,  and  quantization. 

The  first  observation  to  be  made  is  that  the  Doppler  range  measurement  is  subject 
to  group  delay  and  phase  refraction.  Thus  the  time  at  which  the  measurement 
applies,  that  is  the  retarded  measurement  time,  is  equal  to  the  time  of  reception 
minus  one-way  group  delay,  or  simply  the  time  midway  between  pulse  transmission 
and  pulse  reception  (at  the  antenna).  And  the  refraction  term,  which  is  with 
respect  to  the  target  position  at  the  retarded  measurement  time,  represents 
phase,  rather  than  group,  refraction  at  this  position. 

The  bias  error  consists  of  the  combined  instantaneous  values  of  the  range  tracker 
output  error  and  the  serially  uncorrelated  component  of  phase  measurement  error 
at  initialization. 


The  timing  error  has  components  due  to  clock  lag,  differences  in  transmitter 
and  receiver  group  delays,  and  sampler  delay.  Combining  all  terms,  the  syste¬ 
matic  timing  error  has  the  form 


6t  = 


6tc  ♦ 


ATj-  6tR  +  6t$ 


where  6t  denotes  uncompensated  timing  error,  and  the  subscripts  denote  clock 
lag  (C),  receiver  delay  (R),  transmitter  delay  (T),  and  sampler  delay  (S;. 

The  scale  factor  error  consists  of  a  term  due  to  the  speed  of  light  uncertainty 

.Q 

(approximately  10  )  and  a  term  due  to  long  term  clock  frequency  uncertainty 

(approximately  5x20  11  for  a  cesium  clock). 

The  serially  uncorrelated  phase  errors  are  due  to  the  receiver,  clock,  trans¬ 
ponder,  and  sampler.  Ignoring  transponder  phase  noise,  a  typical  worst  case 
phase  error  budget  for  noise,  jitter,  and  quantization  in  CSP  radars  is  presented 
in  Table  5.1.  The  Root  Sum  Squares  (RSS)  error  is  6.2°,  which  is  equivalent 
to  an  error  in  range  of  approximately  0.5  millimeter  for  f  =  5  Gigahertz. 

5.2  Transponder  And  Antenna  Consideration 

The  coherent  transponder  and  antenna  equipment  located  on  the  vehicle  of  interest 
are  vital  subsystems  of  CSP  radar  tracking  operations.  Unsatisfactory  perfor¬ 
mance  of  either  subsystem  can  result  in  erroneous  Doppler  phase  or  frequency 
of  the  signal  replies  to  the  radar. 

Insofar  as  the  transponder  is  concerned,  transponders  under  procurement  for 
Minuteman  III  launches  from  WTR  are  required  by  contractual  specification  to 
remain  frequency  coherent  within  0.6  Hertz  smoothed  over  0  1  second,  but  not 
required  to  be  phase  stable.  Transponder  tests  being  conducted  at  SAMTEC  show 
phase  stability  within  0.001  Hertz  under  some  environmental  condi ti  ^s,  but 
also  display  phase  variation  due  to  signal  strength,  temperature,  interogation 
rate  and  frequency,  and  vibration.  Details  of  the  test  methodology  and  results 
will  be  reported  on  completion  of  the  tests. 
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TABLE  5.1  RECEIVER  RANDOM  ERROR  BUDGET 


RECEIVER  NOISE  AT  SNR  =  20  dB 

4.05° 

CLOCK  PHASE  NOISE 

0.12°  (INTRA-PULSE) 

2.1°  (TRANSIT  TIME:  R  =  5000  nmi) 

TIMING  JITTER  (at  =0.1  Msec) 

4.1°  (ft  =  20,000  ft/sec) 

QUANTIZATION  (8  BIT  I/Q  WORDS) 

0.82° 

RSS  OF  ABOVE  ERRORS  =  6.2° 


* 
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A  basic  factor  in  the  satisfactory  performance  of  the  transponder  antenna  is 
the  stability  of  the  phase  delay  with  respect  to  changes  in  the  aspect  angle 
viewed  from  the  radar,  since  this  aspect  angle  will  change  as  the  missile  moves 
along  its  flight  path  in  addition  to  roll  and  pitch  maneuvers.  An  ideal  antenna 
would  be  the  theoretical  isotropic  radiator  located  at  the  center  of  gravity 
of  the  vehicle,  since  it  would  produce  no  change  in  pnase  with  aspect  angle 
variation.  The  "worst  case"  would  be  two  discrete  radiating  elements  mounted 
diametrically  opposite  on  the  missiles'  cylindrical  surface  wth  equal  power 
applied  to  each  radiator.  If  the  radiators  were  about  8  feet  apart  and  the 
roll  aspect  angle  was  midway  between  the  two  radiators,  a  1.26°  change  in  roll 
aspect  could  result  in  180°  of  Doppler  phase  change  and  signal  strength  nulls 
as  deep  as  -20  dB.  The  Range  Vernier  system  will  interpret  this  change  in  phase 
as  an  erroneous  change  in  range  and  the  CSP  system  will  interpret  the  phase 
change  rate  as  an  erroneous  velocity  change. 

The  closest  approximation  of  the  performance  of  the  isotropic  antenna  commercially 
available  appears  to  be  the  segmented  micro  strip  antenna  design.  This  design 
ideally  will  exhibit  no  change  in  phase  with  roll  aspect  angle,  and  90°  phase 
change  with  a  90°  change  in  pitch  aspect.  It  should  be  remembered  that  the 
antenna  is  a  reciprocal  device  and  any  phase  change  induced  in  the  received 
signal  by  the  antenna  will  have  added  a  like  phase  change  when  the  reply  is 
transmitted. 

The  purpose  of  this  subsection  is  to  emphasize  that  satisfactory  performance 
of  both  the  transponder  and  antenna  subsystems  is  necessary  to  obtain  the  high 
potential  accuracy  achievable  with  the  Range  Vernier  system. 
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6.0 


RANGE  VERNIER  APPLICATIONS 


Doppler  ranges  determined  with  the  Range  Vernier  technique  can  provide  very 
precise  (or  low  noise)  measurements  of  ranges.  Very  precise  range  measurements 
can  substantially  improve  observability  of  trends,  refinement  of  error  models, 
and  accuracy  of  the  estimation  process.  In  turn,  this  can  significantly  improve 
WTR  support  of  weapon  system  test  objectives.  Particular  examples  that  show 
substantial  improvements  gained  from  application  of  very  precise  ranges  are 
presented  in  Sections  6.1  and  6.2  that  follow. 

Section  6.1  shows  certain  results  that  have  been  obtained  by  the  XONICS 
Corporation  with  a  technique  called  Phase-Derived-Range  (PDR).  Additional 
results  and  a  description  of  PDR  are  presented  in  Reference  [3].  The  PDR 
technique  can  be  viewed  as  a  modified  version  of  the  Range  Vernier  concept 
that  has  been  adapted  for  applications  to  skin  tracking,  in  contrast  to  tracking 
a  coherent  transponder  with  a  CSP  radar.  In  view  of  the  similarity  of  PDR  and 
Range  Vernier  techniques,  the  results  shown  in  Section  6.1  are  indicative  of 
contributions  to  data  processing  objectives  that  can  be  anticipated  with  Range 
Vernier. 

A  summary  review  of  several  established  estimation  techniques  and  the  con¬ 
sequences  of  using  Range  Vernier  type  ranges  with  these  techniques,  relevant 
to  the  basic  objectives  of  test  mission  support,  is  presented  in  Reference 
[1].  This  reference  also  includes  an  example  of  Range  Vernier  contributions 
to  Guidance  Analysis  that  is  presented  in  Section  6.2. 

6.1  Operational  Results 

The  operational  results  described  herein  were  presented  by  XONICS,  Inc.  in  the 
paper  cited  in  Reference  [3].  The  results  were  obtained  with  the  PDR  technique 
mentioned  above. 

In  lieu  of  operational  results  for  the  Range  Vernier  that  is  being  developed 
at  SAMTEC,  and  because  of  the  similarity  of  the  two  techniques,  it  seems 
appropriate  that  a  preview  of  results  to  be  anticipated  from  future  applications 
of  Range  Vernier  can  be  gained  from  PDR  results.  Accordingly,  one  of  several 
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case  results  presented  in  Reference  [3]  is  repeated  in  Figure  6.1  that  follows; 
other  results  are  shown  in  the  Reference  but  not  repeated  herein.  The  differences 
AR  shown  graphically  in  Figure  6.1  are  comparisons  of  ranges  from  particular 
radars  with  a  Best  Estimate  of  the  ballistic  trajectory,  from  about  700,000 
feet  altitude  to  about  0.30  second  before  impact.  Comparisons  of  AR  obtained 
with  PDR  versus  usual  range  measurements  collected  from  three  radars,  can  be 
made  from  the  Figure.  The  radars  were  operating  at  VHF ,  S-Band,  and  C-Band 
frequency  levels  respectively.  In  all  cases,  PDR  is  decidedly  superior  and 
improved  observability  is  attained. 

An  important  quotation  from  the  Summary  of  Reference  [3]  is 

"...  .  The  improved  precision  of  PDR  when  applied  to  trajectory  analysis 
results  in  (1)  improved  pierce  points,  (2)  observation  of  altitude  varia¬ 
tions,  (3)  improved  ability  to  estimate  high  altitude  drag  and/or  atmos¬ 
pheric  density,  (4)  improved  observation  of  ionospheric  refraction  effects, 
and  (5)  observation  of  signatures  related  to  RV  configuration  which  have 
diagnostic  as  well  as  discrimination  applications.  ..." 

In  other  words,  the  theoretical  implication  that  extremely  precise  range  measure¬ 
ments  will  improve  observability  of  parameters  of  interest  to  test  missions, 
has  indeed  been  confirmed  with  operational  results.  Moreover,  increased  observ¬ 
ability  gained  with  such  range  precision  can  facilitate  identification,  analysis, 
and  evaluation  of  parameters  and  measurements  that  support  test  mission  objectives 

6.2  Application  To  Guidance  Analysis 

One  of  the  more  difficult  problems  in  ballistic  missile  flight  testing  has  been 
the  analysis  of  guidance  system  errors  at  the  component  level.  Theoretical 
investigations  into  this  problem  have  been  conducted  by  several  organizations 
[4],  [5].  One  measure  of  effectiveness  which  has  been  devised  is  the  recovery 
ratio,  defined  as  the  square  root  of  the  ratio  of  a  posteriori  to  a  priori  error 
variance.  Thus  for  a  given  parameter,  a  lower  recovery  ratio  implies  more 
observability  and  better  recovery  through  the  estimation  procedure. 


4^  R  (feet) 
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FIGURE  6.1  UHF,  S-BAND  AND  C-BAND  PDR  AND  TRACK  RANGE  VERSUS  BALLISTIC  BET 


26 


Analysis  of  Minuteman  guidance  error  recovery  ratios  achievable  by  processing 
conventional  SAMTEC  radar  data  has  generally  shown  that  significant  improvement 
in  guidance  analysis  is  achievable  only  if  radar  accuracy  can  be  improved  by 
some  one  to  two  orders  of  magnitude.  Such  dramatic  improvement  in  the  accuracy 
of  conventional  radar  position  and  Doppler  frequency  measurements  does  not  appear 
to  be  feasible  at  this  time.  However,  similar  analysis  of  Minuteman  guidance 
error  recovery  achievable  by  processing  Doppler  range  data  from  SAMTEC  radars 
predicts  substantial  reduction  of  all  recovery  ratios.  Thus,  with  no  inherent 
improvement  in  radar  accuracy,  but  merely  by  processing  Doppler  phase  rather 
than  Doppler  frequency,  theoretical  studies  predict  substantially  improved 
performance.  If  nothing  else,  this  demonstrates  the  tremendous  loss  of  infor¬ 
mation  which  occurs  when  Doppler  measurements  are  processed  in  the  frequency 
domain  rather  than  the  phase  domain. 

The  results  of  two  simulation  cases  for  the  Minuteman  NS-20  guidance  system 
and  SAMTEC  radars  will  now  be  presented.  The  radar  locations,  measurements 
processed,  observation  intervals,  and  assumed  errors  are  shown  in  Table  6.1. 
and  the  corresponding  recovery  ratios  for  the  guidance  system  error  coefficients 
are  presented  in  Table  6.2.  Observe  that  the  recovery  ratio  of  each  guidance 
term  is  considerably  lower  in  Case  2  when  Doppler  range  measurements  are  processed. 
Furthermore,  although  some  coefficients  are  practically  non-observable  in  Case 
1,  all  are  predicted  to  be  quite  observable  in  Case  2. 

6.3  Realtime  Possibilities 

Although  the  preceeding  examples  and  current  interests  are  primarily  concerned 
with  applications  of  Doppler  ranges  to  post-flight  data  reduction  objectives, 
it  is  germane  to  mention  possible  applications  of  Doppler  ranges  to  realtime 
data  reduction  objectives,  especially  accuracy  improvement  of  impact  predic¬ 
tions  that  are  crucial  to  support  of  range  safety  functions  during  a  test 
mi ssion. 
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TABLE  6.1  RADAR  SCHEDULE 
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STANDARD  DEVIATION 

BIAS _ oR  =  80  FT 

NOISE  0R  =  0.01 


TABLE  6.2  GUIDANCE  ANALYSIS  RECOVERY  RATIOS 


GUIDANCE  SYSTEM 

RECOVERY 

RATIOS 

PARAMETER 

CASE  1 

CASE  2 

0.7250 

0.0065 

PLATFORM  MISALIGNMENT 

0.3143 

0.0040 

0.5075 

0.0152 

0.9887 

0.0303 

GYRO  BIAS  (DRIFT  RATE) 

0.8873 

0.0568 

0.9437 

0.0927 

0.7896 

0.0110 

G-DEPENDENT  DYRO  BIAS 

0.7132 

0.0382 

0.8127 

0.0470 

0.9041 

0.0087 

G2-GEPENDENT  DYRO  BIAS 

0.5534 

0.0087 

0.4799 

0.0091 

0.2369 

0.0005 

ACCELEROMETER  BIAS 

0.3029 

0.0009 

0.2637 

0.0008 

0.4892 

0.0013 

ACCELEROMETER  SCALE  FACTOR 

0.5308 

0.0073 

0.5514 

1 

0.0075 
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Whenever  systematic  error  terms  of  Doppler  ranges  have  been  verified,  it  will 
be  possible  to  remove  these  error  terms  and  provide  corrected  Doppler  ranges. 
In  turn,  very  accurate  range  rate  measurements  could  be  derived  from  poly- 
nominal  smoothing  of  the  corrected  range  data.  This  highly  accurate  range  and 
range  rate  data  could  then  be  used  to  substantially  improve  the  accuracies  of 
estimated  trajectories  and  subsequent  estimates  of  predicted  impact  positions. 
Since  it  is  feasible  to  consider  performing  these  calculations  during  realtime 
data  reduction,  future  considerations  of  possible  applications  of  Range 
Vernier  to  realtime  objectives  is  of  practical  importance  to  WTR  test  support. 
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7.0 


FEASIBILITY  TESTING 


Implementation  and  operational  testing  of  prototype  Range  Vernier  systems  are 
a  major  part  of  the  development  process,  much  knowledge  will  be  gained  as  the 
tests  evolve.  Accordingly,  a  brief  discussion  of  feasibility  testing  is  pre¬ 
sented  herein. 

Relevant  to  the  ultimate  goal  of  providing  very  precise  Doppler  range  measure¬ 
ments  with  modelable  systematic  error  terms,  the  essential  objectives  of 
feasibility  testing  are  to 

(1)  demonstrate  the  operational  ability  of  a  CSP  radar  that  is  augmented 
with  an  (I/Q)  measurement  subsystem  to  provide  the  measurements 
required  (Section  4.3)  by  the  Range  Vernier  technique,  and 

(2)  evaluate  the  ability  to  satisfactorily  model  the  systematic  and 
random  components  of  errors  that  are  present  in  the  Doppler  range 
measurements  obtained  from  the  tests. 

Test  data  will  be  collected  by  tracking  objects  of  opportunity.  Analyses  of 
the  test  data  will  be  made  to  assess  performance  of  the  augmented  system  and 
to  evaluate  various  components  of  the  total  Doppler  range  error.  The  total 
error  will  be  determined  by  comparisons  of  the  Doppler  ranges  with  a  best 
estimate  of  trajectory  (BET)  of  the  object  that  is  tracked.  Accuracy  of  the 
reference  BET  is  crucial  to  evaluation  of  the  total  errors  and  subsequent 
evaluation  of  systematic  and  random  (or  "noise")  error  components.  However, 
initial  estimates  of  the  noise  component  can  be  determined  by  smoothing  the 
Doppler  ranges  with  a  low  order  polynomial,  or  by  comparing  the  ranges  with  a 
smooth  (not  necessarily  a  "best")  estimate  of  the  trajectory. 

Experience  at  SAMTEC  has  shown  that  a  near  earth  satellite  equipped  with  a 
coherent  transponder  antenna  configuration  for  Doppler  tracking,  such  as  GEOS-C, 
is  a  practical  object  of  opportunity  for  various  tests.  The  GEOS-C  satellite 
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offers  (a)  signal-to-noise  levels  sufficient  to  the  needs  of  a  CSP  radar,  (b) 
a  smooth  orbit  that  is  amenable  to  determination  of  an  accurate  BET,  (c)  a 
cross-section  that  permits  skin-tracking,  and  (d)  the  opportunity  to  schedule 
tests  whereby  sufficient  time  can  be  allocated  for  set-up  and/or  calibration 
of  the  radar  immediately  prior  to  actual  tracking  of  the  object. 

Ballistic  missiles  such  as  Minuteman  III,  also  provide  objects  of  opportunity 
for  feasibility  tests.  These  tests  are  important  since  they  will  reflect  the 
actual  environment  of  a  class  of  tests  for  which  a  completed  Range  Vernier 
system  is  intended  to  support.  However,  compared  to  satellite  objects,  several 
comments  regarding  Minuteman  III  objects  are  worthy  of  note;  namely  (1)  repeated 
tests  of  the  same  airborne  transponder/antenna  configuration  are  impossible, 

(2)  the  latter  configuration  is  not  as  stable  as  that  of  GEOS-C,  and  (3)  generally 
speaking,  the  BET  of  a  GEOS-C  satellite  can  be  more  accurate  than  a  BET  of  the 
more  dynamic  Minuteman  vehicle.  Regardless,  of  these  comments,  both  GEOS-C 
type  satellites  and  ballistic  missile  vehicles  provide  good  objects  of  opportunity 
for  feasibility  testing. 

A  note  of  practical  interest  is  that  only  a  relatively  modest  augmentation  of 
the  data  collection  (radar)  and  data  reduction  (software)  capabilities  at  SAMTEC 
is  needed  to  conduct  feasibility  testing.  This  is  because  implementation  of 
Range  Vernier  will  make  full  use  of  current  CSP  radar  capabilities  -  augmented 
with  an  (I/Q)  measurement  subsystem  and  recording  equipment.  Also,  the  analysis 
efforts  will  make  full  use  of  extensive  computer  programs  that  are  currently 
available  at  SAMTEC.  Only  a  relatively  small  computer  program  to  convert  the 
additional  recorded  data  to  Doppler  range  measurements  is  required  to  complete 
software  capabilities  needed  for  analysis  of  the  feasibility  test  data. 
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SUMMARY 


The  Range  Vernier  technique  is  conceptually  straight  forward,  and  can  be 
implemented  by  augmenting  existing  CSP  radars  with  Doppler  phase  extraction 
and  recording  equipment.  Alternative  mechanization  techniques  and  requirements 
for  such  equipment  are  presented.  Development  of  a  prototype  Range  Vernier 
system  is  in  progress  at  SAMTEC,  so  a  discussion  of  feasibility  testing  is 
i ncl uded. 

Resolution  of  Doppler  ambiguities  involved  in  Range  Vernier  and  subsequent 
determination  of  Doppler  ranges  obtained  with  the  technique  are  described.  A 
linear  model  of  systematic  and  random  error  components  of  Doppler  ranges  is 
provided,  including  consideration  of  error  sources  from  an  airborne  transponder/ 
antenna  subsystem  of  a  CSP  radar. 

Substantial  improvements  to  test  mission  support  obtained  from  application  of 
very  precise  range  measurements  are  demonstrated  with  operational  and  theoretical 
examples.  The  operational  examples  clearly  show  that  improved  observability, 
error  recovery,  and  measurement  accuracy  have  been  achieved.  The  theoretical 
example  indicates  how  precise  range  measurements  can  also  facillitate  guidance 
system  analysis.  Applications  of  operational  Doppler  ranges  are  anticipated 
to  provide  very  precise  range  measurements  and  the  improvements  outlined  above. 

In  addition  to  the  post  mission  applications  mentioned  above  possible  application 
of  Doppler  ranges  to  realtime  data  reduction  performed  for  range  safety  is 
di scussed. 
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APPENDIX 


An  example  of  a  method  whereby  Doppler  frequency  ambiguities  can  be  rapidly 
resolved  is  described  herein.  The  method,  previously  mentioned  in  Section  3.3, 
processes  delta-ranges  from  two  sources  with  a  linear  estimator  to  establish 
ambiguity  resolution.  The  sources  are  (a)  Doppler  frequencies  of  the  fineline 
tracking  loop  and  (b)  ranges  from  the  range  tracking  subsystem  of  a  CSP  radar. 
The  linear  estimator  can  be  a  batch  least  squares,  or  a  recursive  Kalman  type 
estimator.  The  former  estimator  is  used  in  the  example  described  below.  In 
either  case,  the  method  is  a  feasible  candidate  for  realtime  or  post-mission 
resolution  of  Doppler  frequency  ambiguities. 

To  begin  discussion,  let  t^  (i  =  0,  1,  . . . ,  N)  be  an  increasing  sequence  of 
time  instants  corresponding  to  Doppler  frequency  measurements  of  successive 
signal  replies  at  the  radar  and  assume  that  a  single  fineline  of  the  pulse 
frequency  spectrum  is  being  tracked  from  t.  through  t...  The  error  6F.  of  F. 
consists  of  a  small  random  component  e.  and  an  unknown  multiple  AN  of  the  pulse 
repetition  frequency  f  .  More  specifically, 

6Fi  =  (F.  -  F.)  =  ANfr  +  e.  ,  (A.l) 

where  AN  is  zero,  a  positive,  or  a  negative  interger,  F.  is  the  true  Doppler 

frequency,  and  the  errors  e.  are  serially  uncorrelated,  with  zero  mean  and 

~2  1 
variance  a  . 

Since  F.  is  related  to  the  range  rate  R(t‘. )  corresponding  to  the  retarded  time 
ti  by  the  equality 

2R(tt) 

F.  =  - - - 

A[1  +  R(tl )/C] 


35 


where  C  is  the  speed  of  light  and  \  is  the  wavelength  of  the  reference 
frequency,  integration  provides  the  delta-range  measurement  ARn 

A*n 

where  the  small  term  R(t.j)/C  has  been  omitted  for  simplicity  of  exposition. 

Accordingly,  ARn  is  a  measurement  of  the  true  delta-range  ARn  =  (Rn  -  Rq)  plus 

an  error  6R  . 

n 

The  error  can  be  determined  by  integration  of  (A.l)  above,  that  is, 


\  /"  F.(t)dt 


n  =  (1,  2,  ....  N), 


6Rn  =  |  /  ANfrdt  +  |  /"  e.dt 

t  t 

0  0 

Consequently, 

6^n  =  I  fr(tn  '  lo)AN  +  ^(n)  ’ 


where 


e(n)  = 


§  / 


e^dt 


To  continue  discussion,  let  Rn  and  Rq  be  ranges  obtained  from  the  range 
tracker  with  errors  en  and  eQ  respectively.  Define  the  delta  range 


Consequently  the  error  6R  of  AR  is 

n  n 


6R  =  (R  -  R  )  -  (R  *  R  )  =  e  -  e 

n  n  n/vo  o'  n  o 
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The  comparisons  yn  can  now  be  considered,  where 

=  A*„  -  4“n  ' 

Since,  AR  -  AR  =  6R  -  6R 
n  n  n  n 

and  6Rn  -  6Rfl  =  |  fr(tn  -  tfl)AN  +  e(n)  -  £„  +  £(), 

the  comparison  y.  (i  =1,  2,  ...,  N)  can  be  expressed  as 

Vi  =  Po  +  X^  +  E.  (A.  2) 

where  =  £Q  ;  X.  =  |  ff(t.  -  tQ)  ;  ^  =  AN  ;  and 

E.  is  a  zero  mean,  serially  uncorrelated  random  error  with  variance  a2  2  variance 
of  £  . . 

For  n  >  2,  with  matrix  notation,  the  least  squares  solution  of  (A. 2)  is 


The  integer  q  for  which  the  absolute  value  -  q  is  minimum  will  be  the 
estimate  of  AN  from  this  example. 
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